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A;, matrix element of four-quark operator relevant to its lifetime is rean- 



^ ' alyzed by QCD sum rule. The new ingredients introduced are that (1) 

^^ ■ more vacuum condensates are considered; (2) different quark-hadron dual- 

^r^ • ity is adopted; and (3) the possible deviation from the vacuum saturation 

^«0 , assumption for the four quark condensates is considered. With ki = 4 the 

Qv I related hadronic parameter r is calculated to be (3.6 it 0.9), and the lifetime 

0^ ' ratio T{Ah)/T{B^) = (0.83 it 0.04), which is consistent with the experimental 



data. 
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I. INTRODUCTION 

Heavy baryon lifetimes provide testing ground for the standard model, especially for QCD 
in some aspects, because they can be calculated systematically by heavy quark expansion 
lU . Theoretically, if we do not assume the failure of the local duality assumption, the heavy 
hadron lifetime differences appear, at most, at the order of ^/toq 0. Recent experimental 
result on the lifetime ratio of the A;, baryon and B meson is [§] 

^ ^ -0.79 ±0.06, (1) 



r 



(50) 



which has shown some deviation from the theoretical expectation. This is one of issues 
concerned with the tests of the heavy quark expansions and has drawn a lot of theoretical 
attention 0-§l, and may imply the potential importance of the 0{l/m\) effect for the above 
heavy baryon and heavy meson lifetime difference P,|9|,p^ 



It is at the order 1/mf that appear four-quark operators, whose contributions to decay 
widths are enhanced by a phase-space factor of (9(167?^) with respect to the leading terms 
in the operator product expansion. Consequently one should include them in the prediction 
for non-leptonic decay rates [|lll]. To the order of l/«^f, the lifetime ratio was calculated as 
follows [|, 

r{K) . , f^lJA,) - ^^IJB) /4CB) 

r(fiO) ^ 2ml < 

HiPiBiirrih) +P2B2{mh) + pseiirrib) + P4<^2ifnb) + [p^ + pQB{mb)]r{mh)} , (2) 

where the term proportional to ^ = (/B/200MeV)^ arises from the contributions of four- 
quark operators. The numerical values of the coefficients cg and j9j's (i = 1 — 6) have 
been calculated in Ref. |Q. fil^Hf,) and /^^(-S) ^ire the averages of 6-quark kinetic and 
chromomagnetic energy, respectively. To the order of 1/ml, the lifetime ratio is 0.98. At 
the scale nib, the values of p^'s are pi = —0.003, p2 = 0.004, p^ = —0.173, p^ = —0.195, 
Ps = —0.012, pq = —0.021. Bi, B2, ci, €2, r and B are the parameterization of the hadronic 
matrix elements of the following four-quark operators, 

(5167^(1 - ^,)qqY{l - l,)h\B) ^ BJlml , 
{B\b{l - 75)gg(l + l,)b\B) ^ B2flml , 
{B\hf,il - l5)taqqY{l - l5)tab\B) = ei/|m| , 

{B\b{l - 75)tagg(l + l5)tab\B) = €2/^3 , (3) 



and 



^ -{Ab\h,il-^,)qqYil-l,)b\Ab)- ^^'"''^ 



2mA, ' ^' ' ' ' ' ' 12 ' 

^ ■(A.IKl - 75)gg(l + l,)b\A,) ^ -B%^r . (4) 



2mA/ ' ' ' " ' ' ' 24 



In Eq. d^), the energy scale for the parameters is at mf,. In Eq. (|^), generally the renormal- 
ization scale is arbitrary, and the parameters depend on it. It can be taken naturally at the 
low hadronic scale to apply the heavy quark expansion. 

These parameters should be calculated by some nonperturbative QCD method. The 



QCD sum rule |T^, which is regarded as a nonperturbative method rooted in QCD itself, 
has been used successfully to calculate the properties of various hadrons. In Refs. |^ and 
0, the mesonic parameters Bi and e, have been calculated by the QCD sum rules within 
the framework of heavy quark effective theory (HQET) [0. The baryonic parameters r 
and B have been calculated in Ref. [[7|. As a result, to the order of 1/m^, the theoretical 
calculation still cannot explain the experiment result (|1|). 

II. THEORETICAL CONSIDERATION 

In this paper, we try to match the experimental data by considering more subtleties in 
the theoretical analysis. As analyzed in Ref. [l^, the lifetime ratio (|l]) depends crucially 



on the value of the baryonic parameter r. Therefore, we shall carry out a reanalysis for 
baryonic parameter r by the QCD sum rule. We note that in the analysis of Ref. [^ , more 
condensates could be included in. They are gluon condensate and six-quark condensate. It 
may be important because the dimension of the gluon condensate is four. In spite of being 
dimension nine, the six-quark condensate may also have significant contribution, because 
the diagram is not suppressed by any loop factor. Experience from the heavy baryon masses 
tells us that the pure tree diagram is important [pISH T^ . Therefore in our calculation, these 
two kinds of condensates will be considered, in addition to those considered in Ref. ^ . 

The conclusion of i? = 1 does not change. It follows from the valence quark approxima- 
tion which is used both in Ref. ^ and in our QCD sum rule analysis. 

To calculate r, the following three-point Green's function is constructed, 

U{uj,uj') = z' fdxdye'''-^-''-y{0\Tf{x)d{0)J"{y)\0) , (5) 



where u = v ■ k and u' = v ■ k'. The Aq baryonic current j" was given in Refs. [|T^-[T8|, 

~f = e^"^ql^C^,ia + b/j)Tq',hl, (6) 

where a and b are certain constants which will be discussed later, /i^ is the heavy quark field 
in the HQET with velocity v, C is the charge conjugate matrix, r is the flavor matrix for 
Aq, 

' ' ^ ^ \ (7) 




In Eq. (^, O denotes the four-quark operator 

O = Kl.^^Kqr^^q . (8) 

Note < At|0|A6 >= — < A^IOIAfo > in the valence quark approximation [0, where 



O = Kl.^Y^qqr^Y^K . (9) 

In terms of the hadronic expression, the parameter r appears in the ground state contribution 
of Il{uj,uj'), 

where the "higher states" denotes the contribution of resonances and continuum. {Aq\0\Aq) 
has appeared in Eq. (^, A = m^^ — uiq and the quantity /a is defined as 

(0|/|Aq) = fAU , (11) 

with u being the unit spinor in the HQET. The QCD sum rule calculations for f\ were 
given in Refs. [|T^-TB]. On the other hand, this Green's function II{uj,uj') can be calculated 



in terms of quark and gluon language with vacuum condensates. This establishes the sum 

rule. For the resonance part of Eq. ([T0|) , we adopt the assumption of quark-hadron duality. 

The calculation of Il{uj,u') in HQET is straightforward. The fixed point gauge |]19| 



IS 



used. The tadpole diagrams in which the light quark lines from the four-quark vertex are 
contracted have been subtracted. The following values of the condensates are used, 

(qq) ~ -(0.23 GeV)^ , 
{asGG) - (0.075 ± 0.015) GeV^ , 
{gqa^.G^'q) = mliqq) , ml ^ 0.8 GeV' . (12) 

Note we have adopted the new value for the gluon condensate from the recent analysis of 



the heavy quarkonium spectrum [20|. However, our final results are not very sensitive to 
the gluon condensate since its contribution is only about 6% of the whole sum rule as shall 
see later. 

Except for the quark gluon mixed condensate, our calculation would be consistent with 
that of Ref. [0] if the gluon and six-quark condensates were omitted. While the calculation 
can be justified if (—a;) and (—a;') are large, however the hadron ground state property should 
be obtained at small (— w) and {—w'). These contradictory requirements are achieved by 
introducing double Borel transformation for uo and a;'. It is defined as 



n — i> oo m ^ OO 

f = ^ fixed f' = ^^ fixed 

n m 




B= lim hm ^-^hrT ^^T^ TT; • (13) 



There are two Borel parameters r and f'. They appear symmetrically, so r = r' = 2T is 
taken in the following analysis. The reason for the factor 2 is similar to that explained in 



the Ref. fET|. 



III. DUALITY ASSUMPTION 



Generally the duality is to simulate the higher states by the whole quark and gluon 
contribution above some threshold energy ujc- The whole contribution of the three-point 
correlator II{(jj,uj') can be expressed by the dispersion relation, 



n(cj,cu') = - 

TT Jo 



/■°° , , Imllfz/, u) 

du / du - 

JO (z/ — c<j)(z/' — a;') 



(14) 



With the redefinition of the integral variables 

Z/ + Z/' 



'^+ 



2 
V — v' 



the integration becomes 



/■OO /"OO /"OO 

/ du l du' ... = 2 du^ 
Jo Jo JO 



'^+ 



du . 



1^+ 



10 Jo 

It is in z/_|. that the quark-hadron duality is assumed [|13 

2 



higher states 



TC Ju), 



du. 



du_ 



Imnfz/, u') 



(15) 



(16) 



(17) 



{u - uj){u' - io') ■ 
This kind of assumption was suggested in calculating the Isgur-Wise function in Ref. |2^ and 



was argued for in Ref. |£2|. As pointed out in p2| , p!3| , in calculating three-point functions 
the duality is valid after integrating the spectral density over the "off-diagonal" variable 
u_ = l(u-u'). 

The sum rule for (Aq|0|Aq) after the integration with the variable u_ is obtained as 



^/Ie--<p(4)<A«l'5|A«) 



duexp 



3(a2 + b^) 



T 



^' + ^'j,8_ «^ „5 



+ Kl 



20487r6 
17(a2 + 62 



^""{giG') + 



8407r6 
5ab 



4871' 



967r2 



{qq) u }- K2 



67r4 

ml{qq)u^ 
ah 



u {qq) 



144 



(18) 



where Ki, K2 are the scale parameters used to indicate the possible deviation from the 
factorization assumption for the four- and six-quark condensates. ki^2 = 1 corresponds to 
the vacuum saturation approximation. ki = (3 ~ 8) is often introduced in order to include 
the nonfactorizable contribution and to fit the data [^]. There is no discussion of K2 in 
literature so we use ^2 = 1. 

As argued in [0, the choice a = b = 1 tends to yield the optimal interpolating current 
for the Ah baryon in HQET. We shall adopt a = 6 = 1 in our numerical analysis. 

The parameters /a, A etc were obtained by the HQET sum rule analysis of two-point 
correlator [p!5HT8[|, 



2 Jo 207r^ TT^ 16 



We have not included as corrections in Eq. (PJ|), because they are also neglected in the 
sum rule for r (ITsD . The values of the parameters are f\ = (2.9 ± 0.5) x 10~^ GeV^, 
A = (0.9 ± 0.1)GeV with the threshold cUc to be (1.2 ± 0.2)GeV and the Borel parameter T 
in the window (0.2 - 0.4) GeV. 

In order to minimize the dependence of the parameters we divide Eq. (|l^) by Eq. ([T9|) to 
extract {Aq\0\Aq) . The variation of the matrix element with Uc for two values of ki = 1, 4 is 
given in Fig. 1. The value oi Uc is (1.2±0.1) GeV. The sum rule window is T = (0.15 — 0.35) 
GeV, which is almost the same as that in the two-point correlator sum rule. We obtain for 
Ki = 4 

(Aq|0|Aq) = (1.6 ± 0.4) X lO-^GeV^^ . (20) 

By taking /^ = 200 MeV, 

r = (3.6 ±0.9). (21) 

If we use fti = 1, we get 

(Aq|0|Aq) = (5.5 ± 1.0) X 10-=^GeV=^ , (22) 

r = (1.3 ±0.3). (23) 

In our numerical calculation, the contribution of the perturbative term is smaller than 
those of the condensates due to small T in the duality window and high dimensionality of the 
spectral function ( z/^ in eq. (18) ). There are several known examples where the perturbative 



part has small contribution in the sum rule ||T6| , p!7|j2^ , without affecting the reliability of 



the results. Nevertheless, a hierarchical structure among various condensate terms exists in 
those cases and the present case. In the present case, because of the cancellation between the 
contributions of the quark condensate and the mixed condensate, as well as the smallness 
of the contributions of the gluon and the six-quark condensates, the influence of the four- 
quark condensate is significant and is enhanced by assuming the deviation from the vacuum 
saturation. 

In Ref. the authors used the following duality assumption, 

1 1 ^ f°° , [°° 1 I Imn(z/, I/') , , 

higher states = — / du I du . (24) 

In their calculation the matrix element (Aq|0|Aq) increases from 0.4 x 10~^GeV^ to 1.2 x 
10~^GeV^ in the working region of the Borel parameter when Uc varies from 1.1 GeV to 
1.3 GeV. In other words, their analysis is very sensitive to the continuum threshold, which 
might imply the above duality assumption is not good. In contrast, the dependence on Uc 
is not so strong in our approach. 

The l/rrife corrections to the above results can be analyzed in principle. While having 
little influence on our above calculation, they formally belong to the 0{l/ml) effects to 



the Ab lifetime. It should be noted that the value of r we have obtained above is at some 
hadronic scale, other than the scale nif,, because we have been working in the framework of 
the HQET, in which the natural scale is ^had "^ f^b- The renormalization group evolution 
of the relevant operators was calculated in Ref . [^,0 • Information on parameters B and r 



at scale m^ is necessary to obtain the lifetime ratio of Eq. (|^). By choosing as{nhad) = 0.5 
(corresponding to fihad ~ 0.67 GeV), Ref. [D gives 



rK)~[1.54 + 0.185(/i;,,rf)]r(/i 



had) 



B{m,) ~ ^^^ . (25) 

1.54 + 0.18E(/i/,,rf) 



That is 



S(mfe) ~ 0.58 , (26) 

r{mi,) ~ (6.2 ± 1.6) (27) 
for Ki = 4 or 

r(mb) ~ (2.3 ±0.6) (28) 



for Ki = 1 from Eqs. ( PT| ) and 

The A{, and B^ lifetime ratio given in Eq. (0) is expressed specifically as 



t(SO) 
for Ki = 4 or 



0.98 - 0.17ei(mb) + 0.20e2(m6) - (0.013 + 0.0225(mb))r(mb) 

(0.83 ±0.04) (29) 



r(Afe 
r(50 



(0.93 ±0.02) (30) 



for Ki = 1. Where the values ei{mi,) = —0.08 and e2{mb) = —0.01 have been taken from 
the QCD sum rules 0. From Eq. (pO]), we see that with the vacuum saturation (ki = 1), 
although r is enhanced by about six times compared to that in Ref. 0, it is still not large 
enough to account for the data Eq. (1). The life time ratio between A^ and B mesons can 
be explained if we also take into account the nonfactorizable contribution of the four-quark 
condensate, as can be seen from Eq. 



IV. SUMMARY AND DISCUSSION 

In summary, we have reanalyzed the QCD sum rule for the A^ matrix element of the four- 
quark operator relevant to the lifetime of A^. Compared to the previous analysis |]^, the new 
ingredients we have introduced are that (i) more condensates are considered; (ii) different 
quark-hadron duality is adopted and (iii) the possible deviation from the vacuum saturation 



assumption for the four quark condensates is considered. Of these ingredients, the last two 
points have more important influence on numerical results than the first point by noting 
the effect due to more condensates is small (~ 10%). The second point is more essential 
than the third point. Note that due to the second point, r is about six times enhanced by 
comparing the result of Eq. (53) with that in Ref. M. And due to the third point, r is about 



three times enhanced by comparing the results of Eqs. (|23| ) and (|2T|) . With fixed duality 
assumption for the sum rule, our result shows that the baryonic parameter r is significantly 
dependent on the nonfactorizable effect of the four-quark condensate and r{mh) = (2.3 ±0.6) 
for Ki = 1 and r{m,f,) = (6.2 ± 1.6) for ki = 4. With this latter value the Af, baryon and B^ 
meson lifetime ratio has been calculated to be t{Ai,)/t{B^) = (0.83 ± 0.04), which is close 
to the experimental data r(A5)/r(i?°) = (0.79 ± 0.06). From the above results we draw the 
conclusion that under the local duality assumption it is possible to resolve the issue of the 
Afo baryon and 5° meson lifetime ratio in the framework of heavy quark expansion. 
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Figure captions 

Fig. 1. Sum rule for (Aq|0|Aq), where O has been given in Eq. (j^). From top to 
bottom these curves correspond to ujc = 1.3, 1.2, 1.1 GeV respectively. The sum rule window 
isT= (0.15-0.35) GeV. 
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